Early fragmentation in the adaptive voter model on directed networks 
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We consider voter dynamics on a directed adaptive network with fixed out-degree distribution. A transition 
between an active piiase and a fragmented phase is observed. This transition is similar to the undirected case if 
the networks are sufficiently dense and have a narrow out-degree distribution. However, if a significant number 
of nodes with low out-degree is present, then fragmentation can occur even far below the estimated critical point 
due to the formation of self-stabilizing structures that nucleate fragmentation. This process may be relevant for 
fragmentation in current political opinion formation processes. 



I. INTRODUCTION 

The defining feature of complex systems is the emergence 
of collective phenomena from the interaction of many parts 
[1]. A vivid example is provided by social (swarm) intelli- 
gence [2]. Crowds of humans, shoals of fish, and even swarms 
of insects are known to solve problems more efficiently than 
any individual on its own [3, 4]. However, only in humans 
social intelligence is also used on a higher level. Among all 
species, only we have evolved the ability to discuss (and de- 
bate) future problems and opportunities, and formulate long- 
term policy. 

Recent media reports have pointed out that a central func- 
tion of the political opinion formation process, the debunking 
of counter-factual opinions, may be starting to fail (see [5] for 
an example from U.S. politics). 

In the political discourse our opinions are exposed to close 
scrutiny and criticism. Well-founded criticism may cause us 
to change counter-factual beliefs and thus promote rational 
decision making. However, because of the stresses involved, 
humans tend to favor discussing with others who share similar 
beliefs. 

With the increasing diversity of offline and online media 
[6, 7] and new media technologies [8], it is becoming easier 
to avoid opposing opinions altogether In particular, the In- 
ternet enables people not only to access but also to publish 
information easily. One of the best examples is perhaps the 
micro-blogging service Twitter, which currently is approach- 
ing 10^ user posts per week. Among this flood of information, 
it is easy to find sources supporting almost every conceivable 
opinion, while avoiding contradicting evidence. 

It thus seems likely that situations develop where a given 
subset of the society (and the media by which it is represented) 
pay attention only to information sources with the same be- 
lief system, thus reinforcing and perpetuating myths that are 
never confronted with opposing views. In this light, one may 
ask whether we are heading for a society that is fractionated 
into groups adhering to internally consistent but mutually ex- 
clusive belief systems. 
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The question when fragmentation in opinion formation pro- 
cesses can occur on a social scale is beyond the scope of clas- 
sical social research, because it focuses on an emergent phe- 
nomenon that may require new conceptual approaches from 
physics. A minimal yet paradigmatic physical model of opin- 
ion formation is the voter model [9, 10]. It describes a network 
of nodes representing agents, and links representing the social 
contacts among them. The agents hold one of two possible 
opinions, which they can change by adopting the opinion of 
their topological neighbors. Due to its similarity to interact- 
ing spin models, the voter model has attracted considerable 
attention in the physics literature and has been studied using 
different interaction geometries, such as regular lattices or het- 
erogeneous networks [11-14]. 

An important extension of the voter model is achieved by 
including homophily, the agents' propensity to discard links to 
opposing neighbors and establish new links to agents holding 
the same opinion [15]. Thus, the interaction network is not 
static but co-evolves with the agent dynamics, as the agents 
rewire their links depending on their opinions. Such networks, 
in which the node and link dynamics co-evolve, are called 
adaptive networks [16-18]. 

The long-term dynamics of the adaptive voter model can 
reach one of two absorbing states: a consensus state, in which 
all agents hold the same opinion, or a fragmented state, in 
which the network breaks into at least two components, which 
are internally in consensus. The adaptive voter model there- 
fore provides a simple framework in which the fragmentation 
of opinion formation processes can be studied [19-23]. In 
this model, fragmentation occurs in a phase transition that has 
been identified as a generic absorbing transition [20]. It can 
be computed analytically using a recently proposed motif ex- 
pansion approach [24] . 

The adaptive voter model has been studied so far on undi- 
rected networks. In the context of opinion formation, how- 
ever, the underlying interactions are often asymmetric. It is 
therefore reasonable to encode "who pays attention to whom" 
as directed links in the interaction network. Directed links 
were considered in previous studies on static networks [25- 
29], but directed adaptive networks have only been investi- 
gated in a generic threshold model for Boolean networks [30]. 

In this paper, we investigate voter dynamics in a directed 
adaptive network, in which both the opinion dynamics and 



2 



the topological change are influenced by the directionality of 
the interactions. The agents can avoid disagreeing with neigh- 
bors by rewiring their outgoing links, whereas they cannot af- 
fect their incoming links. The overall network topology thus 
changes adaptively, while the out-degree of each agent, i.e., 
its number of outgoing links remains unchanged. This en- 
ables us to study the influence of different realistic out-degree 
distributions. 

For sufficiently dense Poissonian out-degree distributions, 
fragmentation occurs if the rewiring of network connections 
exceeds a critical rate, which is consistent with previous re- 
sults on undirected networks [19-23]. However, for scale- 
free out-degree distributions and Poissonian distributions with 
small mean degree, we find that fragmentation can already be 
observed at much lower rewiring rates than in undirected net- 
works. We show that this behavior is due to the nodes of low 
out-degree, which can form self-stabilizing topological struc- 
tures that nucleate fragmentation. 



II. MODEL 

We consider a network of nodes representing agents and 
K directed links representing social contacts. Each node i 
holds a binary opinion (Ji e {A, B}. The direction of links 
indicates the flow of attention between the agents. In other 
words, in our notation we draw links in the direction that one 
would draw the "follows"-links on Twitter. 

We initialize the network as a random directed graph with 
randomly assigned equiprobable strategies. The node states 
and the network topology are then left to evolve according to 
the following rules: In sequential updates, a link i-^j is picked 
at random from the network [31]. If <Ti ~ aj, the link is said 
to be inert and nothing happens. Otherwise, the link is said to 
be active, and it is either rewired (probability p), or an opinion 
update takes place (probability 1 — p). In the former case, the 
node i cuts the link and reconnects to a random node k with 
(Tfc — (Ji. In the latter case, node i switches its opinion cji to 
(jj . We note that the rewiring of links changes only the in- 
degree distribution, whereas the out-degree distribution and 
the average degree (fc) ~ K/N of the network remain fixed. 

In contrast to previous studies of the voter model on static 
directed networks, we do not need to restrict our model to 
networks consisting of a single strongly connected compo- 
nent [25, 28, 32], because the network's component structure 
is affected by the ongoing rewiring of links, which contin- 
uously forms and re-routes paths between different strongly 
connected components. 

Below, we study the proposed model in terms of the density 
n of A-nodes (corresponding to agents holding opinion A) and 
the per-capita densities / and h of active links, / = [A-^B] 
and h = [B— >^A]. Following [20, 21] we characterize the state 
of the network by the magnetization m = 1 — 2n and the 
active link density p = {f + h)/{k). 

In network simulations of the directed adaptive voter model 
one observes qualitatively different types of trajectories: Fkst 
at sufficiently high rewiring rates the network rapidly ap- 
proaches a fragmented state (|m| > 0, p = 0), in which the 
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FIG. 1. (Color online) Typical trajectories from network simula- 
tions. The state of the network is characterized by the density of 
active links p and the magnetization m. The trajectories shown cor- 
respond to networks with Poissonian out-degree distributions with 
(fc) = 4 (red/light gray), (k) = 8 (black), and an out-degree distri- 
bution following Pout (fc) oc (blue/dark gray). The trajectories 
initially drift along a parabola of active states (dotted red, dashed 
black, and dash-dotted blue lines, denoting analytical results from 
Eq. (10)). However, only the black trajectory reaches a consensus 
state, whereas the others eventually collapse to a fragmented state. 
The inset shows a time series of p from the scale-free network shortly 
before fragmentation. The parameters are N = 10" andp = 0.1. 



network breaks into at least two components, which are in- 
ternally in consensus. Second, for sufficiently low rewiring 
rates the network initially approaches a stationary active state 
(|m| <C 1, p > 0), in which the opinions and the topol- 
ogy change continually. Because such active states form a 
parabola in the p-m-plane, the system can drift randomly 
along the parabola until an absorbing consensus state (m = 
±1, p — 0) is reached. These dynamics are closely remi- 
niscent of the adaptive voter model on undirected networks 
[20, 22]. 

In addition to the trajectories described above, the directed 
model can show a third type of behavior not observed in the 
undirected case. Here, the systems drifts along the parabola 
of active states for some time and then collapses slowly to the 
fragmented state (Fig. 1). This can lead to fragmentation sig- 
nificantly below the critical rewiring rate found in undirected 
networks. The delayed fragmentation after the drift along the 
parabola of active states suggests that in the active state the 
network undergoes some slow reorganization that eventually 
leads to the destabilization of the active states. In the fol- 
lowing we investigate the nature of this reorganization and its 
implications for network fragmentation. 



III. ANALYTICAL APPROACH 

The main goal of this paper is exploring the impact of di- 
rectionality of attention on the opinion formation process. For 
this purpose we compare the dynamics of the directed adap- 
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tive voter model to well known results for the undirected adap- 
tive voter model. In the following we refer to these two mod- 
els simply as the directed model and the undirected model, 
respectively. A direct comparison of different models is of- 
ten difficult and may lead to misleading results. We therefore 
compare simulation results of the directed model to two ana- 
lytical approximations that are known to capture the dynamics 
of the undirected model in different limits. In this compar- 
ison an agreement between analytical and numerical results 
indicates that the assumptions made for the undirected model 
are still valid in the directed model, whereas a disagreement 
points to new physics in the directed model that is not ob- 
served in the undirected model. 



A. Moment expansion 

The undirected model was studied extensively by moment 
expansions [20, 22, 33]. Following previous works, we de- 
rive differential equations for the time evolution of so-called 
network moments, namely the densities n, /, and h defined 
above. The change in the density of A-nodes, n, is given 
by the balance between opinion adoption in B-^-A and A-^B 
links, 

n={l-p){h- f). (1) 
The density of A— >B links, /, changes according to 

/ = -pI + (1 - p){[A->A^B] - 2[B^A-^B] 

+ [B^B-^A]-[A-^I^A]-/}, (2) 

where [X— >y— >Z] and [X<— F— >Z] denote the per-capita densi- 
ties of triplets in the network {X, Y, Z e {A, B}). In Eq. (2), 
the first term corresponds to the gain in / due to rewiring, 
whereas the remaining terms correspond to gains and losses 
due to opinion adoption. In an opinion adoption event, a 
node copies a neighbor's opinion via one of its outgoing ac- 
tive links, transforming it into an inert link. This also affects 
all other links connected to the focal node, transforming ac- 
tive links into inert ones and vice versa. The resulting indirect 
change in the density of active links is accounted for by the 
triplet variables. 

The time evolution of the density of B— >A hnks, h, is deter- 
mined by the analogous equation 

h = -ph + (1 - p) { [B->B->A] - 2[Af-B->A] 

+ [A^A-^B] - [B->A-^B] - h}. (3) 

Equations (l)-(3) do not constitute a closed ordinary dif- 
ferential equation (ODE) system, as they involve the triplet 
moments [X^Y^Z] and [X^Y^Z]. In principle, the equa- 
tion system could be complemented by similar equations for 
the triplet moments. These would, however, depend on higher 
moments, such as four-node motifs. An appropriate trunca- 
tion of the expansion is thus necessary in order to obtain a 
closed system of equations. This approximation is referred to 
as moment closure [33, 34]. 



In the following, we express the triplet densities in terms of 
node and link densities using the pair approximation [33-36], 



[X^Y^Z]^[X^Y]{k)^^, 
fixziX^Y^Z] « [Y^X]{qo)^^^, 



(4) 
(5) 



where ^xz = 1 + 5xz accounts for the double-counting of 
symmetric triplets. In the equations, the numbers of triplets 
are approximated as the number of X^Y or Y^X links 
times the average number of Y^Z links connected to a Y- 
node. Here, we assume that the probability of finding an 
outgoing Z-neighbor of a y-node is independent of the pres- 
ence of an X-neighbor of the F-node. Thus the probability 
of finding a given triplet depends on the global link density 
[Y^Z]/{k)[Y]. In (4), each of the (fc) outgoing links of the 
y-node is a Y~^Z link with this probability. In (5), in con- 
trast, the y-node has already been selected by following one 
of its outgoing links. In this case, each of its remaining {qo) 
outgoing links is a Y-^Z link with this probability. The quan- 
tity ((7o) is the mean excess out-degree in the network, which 
can be computed from the out-degree distribution [37]. 

As a further simplification we assume that the mean degree 
of both node types is equal, so that [A— >A] = {k)n — f and 
[B-^B] = (fc) (1 - n) - ft, [20]. We then obtain a closed set of 
ODEs, 



= (h-f), 



(6) 



J^^,mil-n)-h)h_J^ (7) 
1 — n 1 — n n 



h:^-ph + {l~p)\ ^^^^^ '-~h 

[_ 1 — n 

_fA^,mn-f)f _^J^ (8) 
n n 1 — n 

where n = (go) / {k). 

This system of differential equations has a trivial solution. 



h* = r = 0, 



(9) 



corresponding to the absorbing states, in which no active links 
are left. These states can be either fragmented states (0 < 
n* < 1) or consensus states (n* = or 1). 

Additionally, there is a continuum of non-absorbing, active 
stationary states, 

/.*=r=n*(l-n*)(W-(^^^^), (10) 

in which neither consensus nor fragmentation is achieved. 
These states form the parabola in the p-m-plane discussed 
above (Fig. 1). 
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FIG. 2. Transitions of a general active motif (min, rriout, kn, lout)- 
The transition probabilities follow from the link update rule of the 
directed adaptive voter model. New active motifs are created only 
when an opinion update occurs on an incoming active link (second 
transition). In this case the number of incoming and outgoing active 
links (dashed) of the new motif is estimated based on the in-degree 
distribution Pin and out-degree distribution Pout of the underlying 
network. 

In the undirected model the moment expansion is known 
to yield good results when p is far from the fragmentation 
point, but to become less precise for p close to fragmentation. 
To estimate the fragmentation point we therefore resort to a 
motif expansion described below. It is nevertheless instructive 
to consider the critical rewiiing rate 

p, = l-[{k){K+l)]-\ (11) 

which is computed by a linear stability analysis of the ac- 
tive states. For p < pc the active states are stable, whereas 
the absorbing states are stable for p > pc- This estimate of 
the critical point closely resembles the analogous result for 
the adaptive voter model using link update on undirected net- 
works [36]. The main difference is that in the undirected 
case, the parameter k is unknown and therefore usually set 
to unity. This approximation is well justified, because in the 
undirected model, the ongoing rewiring leads to an approxi- 
mately Poissonian degree distribution for which k = 1 [37]. 
In contrast, in the directed model, the outgoing degree distri- 
bution remains fixed, and k has to be considered explicitly. 
In the results presented here, we use values of n that are ex- 
plicitly measured in realizations of the respective out-degree 
distributions. 



B. Motif expansion 

In the undirected model a precise estimate of the transi- 
tion point can be obtained by a motif expansion proposed in 
Ref. [24]. In this expansion we consider an almost fragmented 
network, consisting of two almost isolated components that 
are internally in consensus, but are still connected by a low 
density of "active motifs". In the directed network the ac- 



tive motifs are characterized by their numbers of inert incom- 
ing and outgoing links (min, rrio^t) and active incoming and 
outgoing links (^in,^out)- Following [24], we derive a set of 
balance equations capturing the effect of all possible update 
processes on the densities of active motifs, 

pij^iTi: ^out 1 ^in; ^out) pij^iri: ^out ; ^in; ^out) 

+ 1 PP{'m-in,mout,kn + l,^out) 

I -t + tout 

+ -J — ^"\^_^^, — (1 -p)p{min - i,mout,kn + 1,'out) 
I i ~r tout 

+ ; ^ i PP("^in, m-out - 1, ^in, ^out + 1) 

'in ' ^out ~r ^ 

H ^^^^^ (1 -p)p(^in,^out,TOin,TOout) (12) 

rriin + "T-out 

for mjn > and mout > 1^ and 

p(0, 1, /in, lout) = -P(0, 1, kn, lout) 

+ [1 - p)Pin{ku)Poutilout + I) 

— P(TOin , Tlout , J^in , "■out j 

^^in I "out 

for = and mout — 1- The summation runs over all ac- 
tive motifs (win, TOout, "in, "out) up to maximum in- and out- 
degrees, i.e., over all possible 4-tuples with m,\^ + nin < fcin 
and mout + "out < ^out, where fcin and /cout denote the cut- 
offs. Note that the dimension of the transition matiix grows 
with the cut-off faster than k?' , which has a significant effect 
on the computation time. 

A schematic representation of the transition probabilities is 
shown in Fig. 2. We account for heterogeneous in- and out- 
degree distributions (Pim^out), but assume that the in- and 
out-degree of a node are uncorrected. In the balance equa- 
tions, the fragmented state is obtained as the stationary solu- 
tion containing zero active motifs. The critical rewiring rate pc 
is then extracted from the linear stability analysis of this state 
as the rewiring rate at which the fragmented state becomes 
stable. 



IV. NUMERICAL EXPLORATION OE EARLY 
FRAGMENTATION 

In the following we compare the estimated fragmentation 
points, obtained from the approximations above, with results 
from agent-based simulation of the networks. We first con- 
sider the case of a network with Poissonian out-degree distri- 
bution with (fc) = 8. As a second example we study a network 
with a scale-free out-degree distribution, in which the frag- 
mentation occurs much earlier. We conjecture that this early 
fragmentation occurs due to the presence of a large number of 
nodes with low out-degree, which is then verified in a network 
with Poissonian degree distribution with (/c) = 4. 
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FIG. 3. Fragmentation of a network with Poissonian out-degree dis- 
tribution and (k) — 8. Shown is the absolute value of the magne- 
tization in the final frozen state (top) and the proportion r of sim- 
ulation runs that reaches a fragmented state before fmax ~ 80 000 
(bottom) as a function of the rewiring rate p. Each point is an av- 
erage over 100 simulations. The critical point computed by the mo- 
ment closure approximation (dotted line) overestimates the critical 
rewiring rate, whereas the motif expansion yields a better estimate 
of the true fragmentation point. For the motif-expansion a cut-off of 
kin = feout = 10 was used. For a higher cut-off the estimated critical 
rewiring rate is expected to shift to slightly higher values. 



FIG. 4. Early fragmentation in scale-free networks. The plots are 
analogous to Fig. 3, but describe networks with an out-degree dis- 
tribution Pout(fc) oc Fragmentation occurs far below the esti- 
mated transition points [for the motif expansion a — 2 (dashed line), 
a = 2.25 (dash-double-dotted line), a = 2.5 (dash-dotted line), 
and for the moment closure a — 2 (dotted line)] and extends over a 
wider range. The parameters are A'^ = lO'', fen = fcout = 10, and 
(k) = 5.5665, 2.7042, 1.876 for q = 2, 2.25, and 2.5, respectively. 



analytical approaches capture the dynamics with similar pre- 
cision as in the undirected case. 



A. Poissonian out-degree distribution with (fe) = 8 

We first consider networks with a Poissonian out-degree 
distribution, because this distribution closely matches the dis- 
tribution observed in the undirected model [20]. Starting from 
a random graph with both in- and out-degrees drawn from a 
Poisson distribution with mean (fc), we simulate the full net- 
work dynamics for systems of up to iV = 10^ nodes until a 
frozen state is reached or a maximum simulation time imax is 
exceeded. Time is measured in units of 1/ K, so that K update 
events take place in one simulated time unit. 

The results in Fig. 3 show a relatively sharp fragmentation 
transition at a critical rewiring rate pc ~ 0.79. For p < pc, 
the network reaches a state of global consensus, in which 
all nodes have the same state (|m| — 1). By contrast, for 
p > it separates into two disconnected components of ap- 
proximately the same size, which hold opposing opinions but 
are internally in consensus. These results are strongly remi- 
niscent of the undirected model [20, 22, 24]. 

The analogy between the directed and undirected model 
extends also to the analytical results. As in the undirected 
model the moment expansion overestimates the transition 
point, whereas the motif expansion yields a relatively precise 
estimate. 

The study of the directed model in Poissonian networks 
with (k) — 8 highlights the similarities between the directed 
and undirected networks and provides a basic test for our an- 
alytical approaches. For these networks the directed model 
exhibits the same dynamics as the undirected model and the 



B. Scale-free out-degree distribution 

We now ask how the model behaves for more realistic 
out-degree distributions that cannot be realized in the previ- 
ously studied undirected model. In the following, we consider 
power-law distributions of the form Pout(^) oc A:~", which 
capture the diversity that is observed in a wide variety of so- 
cial applications [38]. For generating networks with power- 
law distributed out-degrees and Poissonian in-degrees, we first 
draw an out-degree sequence of length N from a power-law 
distribution. For each out-degree ki in this sequence, we then 
connect the outgoing links of node i to ki random nodes in 
the network. We explicitly avoid creating nodes without out- 
going links as these would never change their state and act as 
"zealots," trivially preventing the possibility of global consen- 
sus [39]. 

The results in Fig. 4 show that in scale-free networks with 
a e [2,3] the fragmentation occurs much earlier than in the 
Poissonian case. Moreover, the proportion r of networks 
reaching fragmentation now increases gradually with increas- 
ing p. Considering individual simulation runs in detail one 
finds that the networks remain for some time in an active state 
before slowly approaching fragmentation — a behavior not ob- 
served in the undirected model or in the directed networks 
considered in the preceding section. 

The observation that the networks spend some time in 
the active state before fragmenting indicates that these states 
are still feasible at least in the beginning of the simulation 
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FIG. 5. Average in-degree vs. out-degree at the initial (t = 0, 
circles), intermediate {t — 20, squares), and stabilized stage be- 
fore full fragmentation (t — 500, diamonds). Shown is the average 
in-degree of the network nodes vs. their logarithmically binned out- 
degree, averaged over 10 simulation runs. On average, the in-degree 
of nodes with large out-degree decreases with time. The parameters 
are: = 5000 and Pout(fc) oc p = 0.3. 



runs. The mechanism by which fragmentation is reached must 
therefore differ from the mechanism observed in the cases 
studied so far, where fragmentation occurs due to the desta- 
bilization of the parabola of active states in a transcritical bi- 
furcation. 

Notably both the moment and motif expansion seem not 
to capture the different mechanism for fragmentation because 
they significantly overestimate the fragmentation point. In 
the following we call fragmentation well below the estimated 
fragmentation point early fragmentation. The main assump- 
tion used in both approximations is the absence of correla- 
tions between a node's in- and out-degree and between near- 
est neighbors. Their failure thus indicates the appearance of 
coiTelations that are absent or not substantial in the networks 
with Poissonian out-degree distribution. 

Network simulations suggest that early fragmentation is ini- 
tiated by the formation of self-stabilizing structures among the 
agents. To understand the process leading to such structures, 
consider that the networks, even far from fragmentation, are 
partially ordered. On average the number of incoming and 
outgoing neighbors of an agent that share the focal agent's 
opinion will be greater than the number of neighbors that op- 
pose the focal agent's opinion, because the rewiring dynam- 
ics transforms active links into inert ones. This implies that 
if an agent's opinion changes, the agent is likely to experi- 
ence a subsequent loss of incoming links because the majority 
of neighbors now oppose the opinion and rewire their links 
with some probability. In the long run agents that frequently 
change their opinion have lower in-degree than those who 
change their opinion rarely. Therefore, the attention, mea- 
sured in terms of incoming links, focuses on the agents that 
have a low out-degree and thus rarely change their opinion 
(Fig. 5). 

Focusing the attention on agents of low out-degree impedes 



the propagation of opinions across the network. In particular, 
it can lead to the formation of small clusters that have few out- 
going links and hence have a very high resistance to invasion 
of the opposing opinion. In an extreme case small subgraphs 
can form in which all nodes are in consensus and all outgo- 
ing links starting within the subgraph lead to other nodes in 
the same subgraph. Such structures are attracting components 
of the network, which can be entered by following a directed 
path but never left due to the lack of outgoing links. There- 
fore, such subgraphs can never be invaded by the opposing 
opinion. Further, no outgoing links leaving the subgraph can 
be formed because none of the nodes in the subgraph will ever 
rewire an outgoing link. 

We call subgraphs that are hard or impossible to invade self- 
stabilizing structures. The initial formation of such a structure 
through the creation of a state-homogeneous attracting com- 
ponent is a stochastic event that occurs with a small proba- 
bility. However, once such a structure has been formed it 
can grow as other nodes rewire their outgoing links into the 
structure. Nodes of low out-degree can be recruited rapidly 
because only few rewiring events are necessary to rewire all 
of their outgoing links into the self-stabilizing structure. Re- 
cruitment of nodes with more outgoing links takes longer as 
more rewiring events are required. Note that the size of a self- 
stabilizing structure can only increase. In simulations, net- 
works observed shortly before fragmentation are often found 
to consist of two almost disconnected clusters, which are only 
connected by a few nodes of high out-degree. Because of their 
frequent changes of state, these connecting hubs have very 
few or no incoming links. 

For illustration of the mechanism described above an em- 
bedding of a small network at four time points before frag- 
mentation is shown in Fig. 6. The network breaks into two 
almost disconnected clusters. The remaining connections are 
finally formed by a single hub node. The fragmentation has 
been nucleated by the formation of two state-homogeneous at- 
tracting components: a self-referential cycle consisting of five 
nodes of out-degree one and another consisting of two such 
nodes. 

Given the observations above we can explain the shape of 
the trajectories shown in Fig. 1. Because the formation of 
a self-stabilizing structure is a rare event, they are generally 
not present in the initial network. The system therefore ap- 
proaches the parabola of active states, which is in agreement 
with results from the undirected model and the analytical ap- 
proximations for the directed model. However, while the sys- 
tem drifts along the parabola of active states, self-stabilizing 
structures are eventually formed due to the ongoing rewiring. 
As the self-stabilizing structures grow, the permissible range 
for the magnetization shrinks, effectively arresting m as al- 
most all nodes are recruited into the self-stabilizing structures 
(Fig. 7). Because the last nodes to join the structures are "hub" 
nodes with high out-degree, a relatively high density of active 
links can be maintained for some time. Because the hub nodes 
undergo rapid opinion switches, rewiring can only slowly sep- 
arate them from opposing neighbors, which explains the slow 
fragmentation. The switching and rewiring of the hub nodes 
are clearly visible in the time series of the active link den- 




FIG. 6. (Color online) Snapshots of an exemplary network of A'^ = 
100 nodes with out-degree distribution Pout (fc) oc alt = 200, 
250, 350, and 400 (from top to bottom). Two state-homogeneous at- 
tracting components (blue/dark gray) nucleate self-stabilizing struc- 
tures leading to early fragmentation. Note that the nodes with high 
out-degree (represented by the node size) have very low or zero in- 
degree and keep the network connected for a long time. Here, frag- 
mentation is reached at t ~ 413. 



FIG. 7. (Color online) Typical restriction of the permissible mag- 
netization range by self-stabilizing structures. Shown are the upper 
and lower bounds for m (dashed red and dotted blue lines) computed 
from the total size Sa and of the stable A- and B-components, 
i.e. maximally connected components in which each node has out- 
going links to nodes in the same state only. The magnetization is 
confined to the shaded area and effectively arrested eventually at a 
fixed value. The aiTow marks the point when only a single remaining 
hub node connects the two parts of the network (cf. Fig. 6). Here 
fragmentation is reached at t « 524. The jumps occur when the last 
links leaving a self-stabilizing structure are rewired to nodes within 
the same structure, so that it is fully stabilized and can be detected. 
The parameters are A'^ — 500, Pout(fc) oc k~'^, andp — 0.3. 

sity p. In the inset in Fig. 1, this is shown for the last 10^ 
time units before fragmentation. Here, one of the two remain- 
ing connecting hubs detaches from one of the components and 
the final hub still switches several times before eventually also 
separating. 

Summarizing the observations above, we conjecture that 
early fragmentation is initiated by the formation of self- 
stabilizing structures among nodes of low out degree. We em- 
phasize that contrary to most dynamical phenomena observed 
in scale-free networks, the dynamics of interest is generated 
primarily in the nodes of low degree. Nodes of high degree 
still play an important role as they are the last nodes to con- 
nect the separating components and thus determine the time 
of fragmentation. This mechanism is not captured by cur- 
rent analytical approaches, because it relies on the build-up of 
negative coiTelations between the in-degree and out-degree of 
nodes that is neglected in previously proposed approximation 
schemes. 



C. Poissonian out-degree distribution with low (fc) 

Because the mechanism postulated above relies on the for- 
mation of coiTelations, one can perform a simple test by con- 
sidering a system in which these coiTelations are removed by 
an additional rewiring process. However, such a test is diffi- 
cult in scale-free networks for two reasons: First, because of 
the constraints in scale-free topology it is well-known that it 
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to the model used so far, except that when an inert link is cho- 
sen, this link is also rewired to a randomly chosen target node 
that is in the same state as the source. 

The model variant in which the rewiring of inert links is 
switched on shows no evidence for early fragmentation (see 
Fig. 8). Fragmentation occurs in a relatively sharp transition 
at a critical rewiring rate that is consistent with the estimate 
from the motif expansion. We emphasize that the rewiring 
of inert links neither introduces nor destroys active links. It 
therefore has no direct impact on fragmentation. However, 
rewiring inert links prevents the build-up of correlation be- 
tween the in-degree and the out-degree of nodes and thereby 
inhibits the formation of self-stabilizing structures. The ab- 
sence of early fragmentation in a model where these corre- 
lations are removed confirms the causal relationships postu- 
lated above. We therefore conclude that in directed adaptive 
networks the slow build-up of negative correlations between 
in-degree and out-degree can initiate early fragmentation by 
leading to the formation of self-stabilizing structures. 



FIG. 8. Fragmentation of networks with Poissonian out-degree 
distribution and low (k), in which the inert links are also rewired 
(squares) or are not rewired (circles). Early fragmentation is clearly 
visible, although the shown averages over 100 simulation runs are 
still rather noisy due to the highly stochastic nature of early frag- 
mentation. Shown are the absolute value of the magnetization (top), 
the proportion of fragmenting simulation runs (middle), and the cor- 
relation coefficient between the in- and out-degree of the nodes in 
the final state (bottom). The parameters are A'^ = 10*, (k) = 4, and 

^in ~ koxit ~ 10. 



is difficult to remove correlations in scale-free networks com- 
pletely, and second, because of the presence of nodes of very 
high degree, fragmentation takes a long time, making numer- 
ical studies of fragmentation tedious. 

Our reasoning above predicts that early fragmentation 
should be observed also in directed Poissonian networks with 
sufficiently low mean degree. In the present section we there- 
fore consider a Poissonian network with a mean degree of 
4, which avoids the difficulties encountered in scale-free net- 
works. In this section we show (a) that this network exhibits 
early fragmentation and (b) that the early fragmentation can 
be avoided by an additional rewiring mechanism that destroys 
the correlations implicated in the formation of self-stabilizing 
structures. 

Simulation results for the network described above are 
shown in Fig. 8. The figure shows clear evidence of fragmen- 
tation well below the estimated fragmentation point. Further, 
this early fragmentation is accompanied by the build-up of 
negative correlations between the in- and out- degrees of the 
nodes. This confirms our previous observation that attention 
focuses on those nodes who pay little attention to others them- 
selves, as we have also seen in the scale-free case in Fig. 5. 

To verify that the correlation described above is the cause 
and not a symptom of the early fragmentation, we now con- 
sider a different variant of the model. This variant is identical 



V. CONCLUSIONS 

In the present paper, we have investigated an extension of 
the voter model on adaptive networks that takes the direction- 
ality of the interactions among the agents into account. We 
found that our model can transition to a fractionated state for 
rewiring rates that lie much below the critical value estimated 
using analytical approaches known to work well in undirected 
models. We discovered that fragmentation occurs due to a 
novel mechanism that depends inherently on the directed na- 
ture of the links. This early fragmentation occurs when agents 
focus their attention on those who are steady in their opinion 
because they pay attention only to few sources of information. 
In this case self-stabilizing structures can form that nucleate 
fragmentation. 

Our results illustrate that directed networks can exhibit new 
physics not observed in their undirected counterparts. Espe- 
cially in the investigation of opinion formation processes, the 
often directed flow of attention should therefore be taken into 
account in models. 

In the context of real- world opinion formation processes the 
mechanism described here may constitute a threat but also an 
opportunity. Early fragmentation maintains diversity of opin- 
ions, it may thus aid the survival of counter-factual myths, but 
also of legitimate and well-founded views of minorities. It is 
conceivable that this mechanism may be employed in the fu- 
ture to adjust the perceived degree of controversy in controlled 
environments. For instance, in online discussion boards or 
music recommendation systems the underlying software can 
in principle control which posts are displayed to whom. It 
could thus encourage rewiring of attention that either facili- 
tates or inhibits early fragmentation. 

Even the adaptive directed voter model paints a highly sim- 
plified picture of real-world opinion formation processes and 
thus must be considered as a toy model. Therefore, investiga- 
tion of more realistic models is an important goal for the fu- 
ture. Based on the results and analysis presented in the present 
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paper we believe that the mechanism of early fragmentation 
will be observed whenever directed attention is focused pref- 
erentially on agents that change their opinions at less than 
average rate. We therefore expect that early fragmentation 
should be robust to future refinements of the model. 

A key ingredient that is missing in our present model is 
novelty. Here we considered only the exchange of opinions 
regarding a single well defined question, whereas in reality 
many discussions are enriched by the constant inflow of new 



ideas. We have shown that homophily favors connecting to 
poorly informed agents and thereby promotes early fragmen- 
tation, whereas curiosity would favor connecting to well in- 
formed agents and thereby hinder early fragmentation. In this 
light, novelty, whether in the form of true innovation or arbi- 
trarily changing fashions may play an important role in pre- 
venting social fragmentation. 

The authors thank G. Demirel for fruitful discussions. The 
work of CH was partially supported by the National Science 
Foundation under Grant No. PHY-0848755. 
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